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Power Corrections to e + e Dijet Event Shapes 
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Abstract. We discuss a class of event shapes for e + e~ dijet events that include the thrust as a special 
case. Large logarithmic corrections to the corresponding cross sections can be resummed to all logarithmic 
orders at leading power. However, irrespective of the order up to which the perturbative expansion is 
calculated, it has to be supplemented by nonperturbative corrections due to its at best asymptotic nature. 
We find that the leading power corrections are universal for the class of event shapes discussed here. Based 
on these findings, we provide sample numerical predictions for the distributions of the new event shapes. 

PACS. 12.38.Cy Summation of perturbation theory in QCD - 13.87.-a Jets in large-Q 2 scattering 



1 Motivation 

Event shapes are generalizations of jet cross sections that 
describe the distribution of radiation in the final state 
PP. They are thus sensitive to both the underlying short- 
distance hard scattering and to long-time hadronization 
processes. In the following, although not always stated 
explicitly, we will consider event shapes in e + e~ annihila- 
tion, where long-time effects only affect the final state. 

The hard scattering is responsible for the basic form 
of the distribution of radiation. For example, the contri- 
bution of a two-jet event to an event shape cross section is 
distinctly different from that of a three-jet final state. The 
distribution of radiation stemming from the short-distance 
process can be computed within perturbation theory in 
terms of quarks and gluons. Originally, this was one of 
the main motivations to study event shapes, to test QCD 
III 

The correctness of QCD was quickly established by 
confirming the main features of event shape distributions 
PQ. Despite this success, perturbation theory cannot fully 
account for their exact form. This is due to hadroniza- 
tion effects which enter as nonperturbative (NP) power 
corrections. The main, although not only, effect of the 
hadronization process is to widen the distribution of radi- 
ation in the final state, and thus to shift the peak of the 
event shape distribution away from the narrow-jet limit |21 
|2j • For the description of this shift it suffices to consider 
the first power correction, typically proportional to l/Q, 
where Q is the center of mass (cm.) energy. In general, 
this gives an accurate picture of the energy dependence of 
the average values of a variety of e + e~ event shapes with 
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only a small set of parameters that are fit to experiment 
(see 0] and references therein). For differential distribu- 
tions, however, one needs to take into account more than 
the first power correction, which leads to the introduction 
of nonperturbative shape functions 0. Perturbation the- 
ory imposes constraints on these shape functions such that 
good agreement with experiment can be achieved, again 
with only a small set of NP parameters [S| • 

The present state of the art is perturbative computa- 
tions up to next-to-leading order (NLO) in the strong cou- 
pling, with large logarithmic corrections due to soft and/or 
collinear radiation resummed to next-to-leading logarithm 
(NLL). Supplemented with the aforementioned set of non- 
perturbative parameters, many e + e~ event shape cross 
sections can be predicted with impressive accuracy. One of 
the main applications is thus the precise determination of 
the strong coupling 0]. Nevertheless, as we will illustrate 
below, there is still more to be learned from the theoretical 
study of such event shapes, aside from the computation of 
yet higher orders or their numerical evaluation. 

In Ref. 0, the authors, together with T. Kiics, intro- 
duced a class of event shapes r Q , depending on a contin- 
uous, real, parameter a, that includes as special cases the 
familiar thrust 8 and jet broadening 9 . As we will review 
below in Sectional the introduction of this parameter a al- 
lows us to analytically study a whole range of event shapes 
simultaneously, and from its variation deduce a variety of 
consequences for both the perturbatively calculable contri- 
butions [7] , and the long-distance corrections [lUllllj . The 
factorization of long- and short-distance effects results in 
a scaling rule that relates the power corrections within 
the considered class of event shapes, which is described in 
Sec. |31 We provide numerical examples, and conclude by 
summarizing present and possible future perspectives. 



Carola F. Berger, George Sterman: Power Corrections to e + e Dijet Event Shapes 



For sake of brevity, we only quote here the main fea- 
tures of the analytical results obtained in Refs. and 
More details, an account of various technicalities, and fur- 
ther references can be found there. Here we show some 
new numerical results. 



2 Resummation of Large Logarithms 

We consider the following distributions of event shapes r a 
that weight final states N in e + e~ annihilation processes 
in the narrow two-jet limit, 

\M(N)\ 2 S(r a ~r a (N)), (1) 



iV 



with a center of mass energy Q 3> ^Iqcd- We sum over 
all final states N that contribute to the weighted event. 
M(N) denotes the corresponding amplitude. The weight 
functions defined in Ref. [7] for a state iV are 



^ all i£N 



(2) 



where pi± is the transverse momentum of particle i rela- 
tive to the thrust axis, and r}i is the corresponding pseudo- 
rapidity, r\i = In cot (Oi/2). The thrust axis [Hj is the axis 
with respect to which the above expression is minimized 
at a = 0. Parameter a is adjustable, — oo < a < 2, and 
allows us to study various event shapes within the same 
formalism. The case a = in Eq. J3J) is essentially 1 — T, 
with T the thrust while a = 1 is the jet broadening. 

In the two-jet limit, r a — ► 0, and @J has large correc- 
tions in In (1/to) that have been resummed to all logarith- 
mic orders at leading power for a < 1 in |7j. For a ~ 1, 
recoil effects have to be taken into account, as was pointed 
out for the broadening (a = 1) in Ref. |12| . 

In the following we will quote the result of the resum- 
mation of large logarithms of r a in Laplace moment space: 



a (u, Q,a) = / dr a e 
Jo 



-VT a da(T ai Q) 
dr a 



(3) 



Logarithms of \/r a are transformed to logarithms of v. 

We refrain from quoting the slightly more complex for- 
mula for the resummed cross section valid to all logarith- 
mic orders 0, and only give the result at NLL, 



J u J Pi 
L u 2 Q 2 



- B(a 

2 v ' 



(4) 



The resummation is in terms of anomalous dimensions 
A(a s ) and B(a s ), which have finite expansions in the 



running coupling, A(a s ) — Y^=i ^ (a s /f) n , and sim- 
ilarly for B(a s ), with the well-known coefficients A^ = 
C F , BW = -3/2 C F , AW = l/2C f [CA(67/18-7r 2 /6)- 
10/9 TpNf], at NLL accuracy, independent of a, where 
C F = 4/3, C A = 3, T F = 1/2, and N f denotes the num- 
ber of flavors. At a — we reproduce the NLL resummed 
thrust cross section |13) . 



3 Universality of Power Corrections 

As it stands, the perturbative cross section, Eq. (TJ} is ill- 
defined for small values of r a . This is due to the at best 
asymptotic nature of the perturbation series, which man- 
ifests itself in @ as an ambiguity in how to treat the 
singularities in the running coupling. At r a ~ A QCD /Q 
nonperturbative corrections become dominant. Neverthe- 
less, due to the quantum mechanical incoherence of short- 
and long-distance effects, one can separate the perturba- 
tive part from the NP contribution in a well-defined, al- 
though prescription-dependent manner (see and refer- 
ences therein). 

Following Refs. we can deduce the structure of the 
NP corrections by a direct expansion of the integrand in 
the exponent at momentum scales below an infrared fac- 
torization scale k. We rewrite Eq. as the sum of a 
perturbative term, labelled with the subscript PT, where 
all p± > k, and a soft term that contains all NP physics: 



In 



1 



a(u,Q,a) 



= 2 



Q~ 



Pi 



dpi 
Pi 



A{a s {p!_)) 



p±/Q 



du 
u 



B-terni 



= In 



Pl/Q 2 
~ ^pt (v, Q, k, a) 



n — 1 n 



(5) 



We have suppressed terms of order (D(v/Q 2 



7Q 2 



which include the entire J3-term of Eq. (TJJ, as indicated. 
Introducing the shape function as an expansion in powers 
of v /Q, from the expansion of the exponent in the second 
equality of with NP coefficients A n (fc), we arrive at, 



a (v, Q, a) = iTp T (v, Q, k, a) / Q:NP [ — , k ) , (6) 



n=l 



(7) 



We find the simple result that the only dependence on a 
is through an overall factor 1/(1 — a) which leads to the 
scaling rule for the shape function [10111] : 



(8) 



/a.NP ( Qi K ) ~ 



/o,NP [ q, « 
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Fig. 1. Differential distributions (1 / o tot )da / dr a for a = 0, 
and a = —0.5 at Q = 91 GeV. The solid lines are the full 
cross sections (for a = the output of PYTHIA), the dashed 
lines the perturbative contributions at NLL/NLO. The data 
are taken from 1161. 



Given the shape function for the thrust at a = at a spe- 
cific cm. energy Q, one can predict the shape function and 
thus from Eq. © the cross section including all leading 
power corrections for any other value of a. 

Before showing some numerical results, we want to 
comment briefly on the main assumptions that go into 
the above derivation. Our starting point is the NLL re- 
summed cross section (@J, which describes independent 
radiation from the two primary outgoing partons. Corre- 
lations between hemispheres are neglected, although they 
are present in the resummed formula valid to all logarith- 
mic orders [7]. However, it has been found from numerical 
studies that such correlations may indeed be unimportant 
0. Furthermore, we assume in the separation of pertur- 
bative and NP effects ©, that long-distance physics has 
the same properties under boosts as the short-distance ra- 
diation. Success of the scaling (jHJ) would thus indicate that 
NP processes are boost- invariant as well 



4 Numerical Predictions 

In 22 we showed some examples of the shape function in 
moment space, computed from the perturbative expres- 
sion matched to fixed order calculations at NLO with 
EVENT2 23] , and from the output of the event generator 
PYTHIA 23 m the absence of data analysis for values 
of a different from the thrust (a = 0). Here we show re- 
sults in momentum space, after numerical inversion of the 
Laplace transform. 

In Fig. ^ we show our prediction for the differential 
distribution for a = —0.5 in momentum space. This plot is 
obtained from Eq. (jSJ, with the shape function determined 
via (JSJ from the one for the thrust, found in 2U- For 
further technical details we refer to 23- We also show 
the differential distributions computed from perturbation 
theory, and the output of PYTHIA for a = 0. As shown, 



PYTHIA fits the data for the thrust well, thus we take its 
output to determine the shape function for a = 0, instead 
of fitting a function to the data. 



5 Summary and Outlook 

We have illustrated how one can test specific properties of 
long-distance physics with the help of perturbation the- 
ory, if supplemented by experimental information. Suc- 
cess or failure of the scaling will provide information 
about the importance of interjet correlations and about 
the boost properties of long-range interactions in e + e~ 
dijet events. This additional information about nonper- 
turbative physics in turn could, for example, be helpful in 
precision measurements of a s . We hope that experimen- 
tal tests of our predictions will be carried out in the near 
future. 
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